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Abstract 

While the usual harmonic oscillator potential gives discrete energies in the non- 
relativistic case, it does not however give genuine bound states in the relativistic case if 
the potential is treated in the usual way. In the present article, we have obtained the 
eigenfunctions of the Dirac oscillator in two spatial dimensions, adapting the prescription 
of Moshinsky. 
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The study of quantum harmonic oscillators has received considerable attention and 
hence it is of intrinsic interest to extend the model to the relativistic domain. Moshinsky 1 
formulated a novel prescription of introducing the interaction in the Dirac equation, 
which, besides the momentum is also linear in co-ordinates and subsequently showed 
that the non-relativistic form of the interaction reduces to one of the harmonic oscillator 
type. While in literatuer such a system is called the 'Dirac Oscillator', we suggest the 
name "Moshinsky-Szczepaniak Oscillator", after its proponents. 

Several authors have addressed the Dirac oscillator in one space dimension. Titchmarsh 2 
analysed the relativistic harmonic oscillator problem using the Green's function tech- 
nique. Nogami et.al. 3 have pointed out the interesting differences in the coherent states 
of the Dirac Oscillator and the non-relativistic harmonic oscillator. The relativistic ex- 
tension of the one-dimensional oscillator in the Lagrangian formalism omegaas developed 
by Moreau 4 and others. While Moshinsky 5 has analysed the one-body relativistic oscilla- 
tor using Group Theory, Villalba 6 has dealt with the angular momentum operator in the 
Dirac equation. Dominguez-Adame 7 has anaylsed the one-dimensional Dirac Oscillator 
with a scalar interaction and shown the absence of Klein paradox. 

An interesting framework for discussing the Dirac Oscillator is a 2+1 space-time. 
Presently we explore the two-dimensional Dirac Oscillator and expose some special fea- 
tures not displayed by one-dimensional systems. 

The Dirac equation in two dimensions 8 for a free particle would read 

Eip = (c (a x p x + a y p y ) + /3m c 2 )^ , (1) 

with m denoting the rest mass of the particle and a x , a y and (3 representing the standard 
Dirac matrices. 

The prescription of Moshinsky is extended to two dimensions and the interaction is 
introduced as follows 

Px^Px- if3m wx (2) 
Py^Py- i(3m wy (3) 

The two dimensional Dirac equation may now be written as 

d . „ „ ... 9 



(3E+ihc(3a x — — +icBa x 8mowx+ihc8a v — +ic(3a v (3mowy —8 2 m§(? 
ox ay 



^ = (4) 



It is convenient to introduce the following representation in terms of the Pauli matrices. 

1 J)' a v = °v = (°i ~o) and/^ = a,= (J ^ 

Writing ip(x) = ( ^ V the matrix form of Eq. (4) is 
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/ E — m c 2 

^4hc-^ + foc-jj- — icm wx + m cwy 



ihc^ + hcJ^ — icm wx—m cwy\ (tpi 



'dy 

—E — m c 2 



J 



'0\ 



(5) 



The spinor equation may be written as a system of two first order coupled differential 
equations 



(d d \ 

ihc— + hc— — icm wx — m cwy \ ip 2 — 
ox ay ) 



(6) 



(d d \ 
—iUc— + hc— icniQWX+moCwy ) ipi — [E + m c 2 ) ip 2 = 
ox ay J 

Equation for ipi: 

The equation for ipi is obtained by using Eq. (7) in Eq. (6). 
On simplification we obtain 



(E 2 - m y)^ 1 + {n 2 c 2 ^- 2+ n 2 c 2 92 



2 2 2 2 

m c w x 



2 2 2 2 

m c w y 



dx 2 2 dy- 

+2nc m w + 2m c wL z > ipi — 



where L z = xp y — yp z is implied. 



(7) 



(8) 



d d 

Using p x = and p y = —ih-Q-, it is straightforward to check that the above 

equation may be written in an elegant form as 



hw 



E 2 -m 2 c^ 
2moc 2 



(9) 



Apparently, the first two terms are those that appear in the Hamiltonian of a non- 
relativistic 2D harmonic oscillator. This justifies why the 'potential' is called the Rela- 
tivistic Oscillator Potential. The fact that Dirac particles remain bound by this inter- 
action suggests the absence of Klein paradox. 

Equation (9) may also be written as 



d 2 



+ 



d 2 



rriniu 2 



v dx 2 dy 2 



2m n 



hw 



4>i 



2m 



wL E 2 - m 2 c 4 
2m,QC 2 



V'i (io) 



It is convenient to use plane polar co-ordinates (p, 0) to obtain the solution of this 
equation. 
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Writing p 2 = x 2 + y 2 , and expressing y 2 as 

1 d ( d \ Id 2 



V 2 = -77- I Pw: I + 



pdp \ dp J p 2 d(j) 2 

the wave function may be expressed as 

= i2 (p) $ (0) = R (p) e 1 ^ (11) 

where m = 0, ±1, ±2, ±3, .... is the angular momentum quantum number. 
Equation (10) takes the form 

1 d ( dR\ m 2 ^ . , ~mlw 2 p 2 ^, . 2m hw ^ , . 



2m . _ , . 2m /.E 2 - mlc 4 



mw ^(p) + » £_ i?(p) = (12) 



Multiplying by p 2 , the above equation becomes 



.2" " . p ^ + (Vp 2 - m 2 ) i? - ^p 4 i? = 0. (13) 



<ip 2 dp ^ 1 h 



Here, 



2 2m 



fc 2 = 



ft 2 



( E 2 — m 2 c 4 " 



(14) 



It is easy to check that k has the dimensions of inverse length and thus kp is a dimen- 
sionless parameter. 

In what follows we use the notation z = — f~~P 2 - m terms of the variable z, the above 
equation may be written as 



,d 2 R dR 1 



where 



dz 2 

fci = A; 2 



+ 2— + - f fcjz - m 2 - z 2 } R (z) ■ 
d? -A L J 

(£ 2 - m 2 c 4 ) 



ft 



m w 



2(m + 1) + 



m c 2 hw 



(15) 



(16) 



While the first term refers to the oscillator part, the second terms refers to the kinetic 
energy of the particle. We try solutions of the form 



, , — z m . . 

R(z) = e * z * q>i_(z) . 



(17) 
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The unknown function <pi (z) should be a constant for z — >■ and should guarantee 
normalisation. Equation (15) may be written in the standard form as 



d 2 

Z 



_ + (6- z )_-^ 1 = (18) 



with b = m + 1 and a = ^(ra + 1 — We identify this differential equation as 
the Kummer's equation 9 whose solution can be expressed in terms of the confluent 
hypergeometeric functions. The only admissible solution is M (a, b, z). The other linearly 
independent solution U (a,b,z) w is rejected since it is irregular at infinity. 

We now write the solution as 

=e— p ^p 2 j M(-(m+l-^),m+l,^-p 2 j (19) 
The complete wavefunction is written as 



* = ^-e"^* ((^) >)" M (i (m + 1 - 4- 1, Of/ 



(20) 



where A is the normalisation constant. Further, as is well-known, when a = —n, the 
hypergeometric series terminates and defines a finite polynomial of n th degree. Hence 



l \ m 



^ = ^e»V^V((^)%) M(-„-l, m+ l,^) (21) 

The non-negative integer n, which arises from the boundary condition that the wave 
function vanishes as p — > oo, also quantizes the energy eigenvalues. 

Equation for ip 2 

The equation for the small component of the Dirac wave function is obtained in an 
analogous manner. Eliminating ipi in Eq. (7) first, and then going through similar steps 
as before, we obtain 



.2 x , E 2 -m 2 c i 



Px + Py 



2m 2m 



+ -ttiqW 2 (x 2 + y 2 ^j + hw 



1p 2 



wL z + 



2m c 2 



,{j 2 (22) 



which resembles Eq. (9). Adapting the same procedure as before it is seen that the 
eigenfunctions can be expressed as 

(1 \ rn 
{^f)%) M(-„,m + l,I^V). (23) 
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The Energy Spectrum 

The quantisation of energy demands the vanishing of the eigenfunctions at infinity, 
which comes from the fact that the first argument of M(a, b, z) be equal to a negative 
integer or zero 8 . We obtain from Eq. (20) 

= -n - 1 (24) 
or equivalently 

E 2 - m 2 c 4 = 4 (n + 1) m c 2 kw. (25) 

As is expected for central potentials, the energy eigenvalues are independent of the 
quantum number m, as a consequence of rotational symmetry. 

Non-relativistic limit 

The non-relativistic limit is obtained by setting E = ttiqC 2 + e nr and considering 
e nr « uiqc 2 . It is straightforward to check that Eq. (25) may be written as 



E = m,Q(? 



( 4 (n + 1) hw 



m c 2 



(26) 



Taylor expansion of the above would give us 



9 / x + 2 (n + l) 2 h 2 w 2 
E ^m c 2 + 2{n + l)hw ^ '- . 27 

TTIqC 2 

It is thus seen that the first term corresponds to the rest energy of the particle, the 
second term refers to the non relativistic harmonic oscillator energy and the third is the 
relativistic correction term. 

Results and Discussion 

The harmonic oscillator problem in quantum mechanics has far-reaching consequences. 
While the non-relativistic harmonic oscillator problem is addressed by solving the Schrodinger 
equation directly or by the well-known matrix formulation method or the operator 
method, the relativistic harmonic oscillator requires an altogether different formalism. 
The prescription of Moshinsky which describes the one-dimensional Dirac oscillator may 
well be extended to two space dimension. It is seen that the eigenfunctions are expressed 
in terms of regular confluent hypergeometric functions, special cases of which are the 
well known Hermite polynomials. More importantly the eigenenergies have the appro- 
priate non-relativistic limit. The Dirac oscillator is thus a relativistic generalisation of 
the quantum harmonic oscillator. 
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Unlike the non-relativistic oscillator, where the energy levels are discrete and equi- 
spaced, the relativistic oscillator no doubt has discrete energies, but unevenly spaced 
levels. Dirac oscillator has potential applications in models of quark confinement in 
particle physics. 
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